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I.  INTRODUCTION 


In  this  report,  on  the  basis  of  elastic  dynamic  analyse,  analytical 
means  have  been  developed  to  design  a  fuselage  for  a  given  blast  intensity 
or  to  determine  the  maximum  blast  a  given  fuselage  can  withstand.  Charts 
in  nondimensional  form  have  been  developed  determining  the  critical  blast 
intensity  for  a  given  skin  size.  From  the  skin's  response,  the  loading  on 
the  stringers  has  been  obtained  and  is  given  graphically  for  various  values 
of  a  nondimensional  parameter  for  peak  pressure  intensities  of  PQ  =  2  psi 

'and  PQ  =  1  psi,  and  their  utilization  is  presented  for  different  skin  sizes. 

v 

The  response  of  the  stringers  was  formulated^ and  the  maximum  bending 
moments  produced  in  a  typical  stringer  for  peak  pressure  intensity  of 
Po  -  2  psi  and  PQ  =  1  psi,  for  different  panel  sizes,  are  presented  in 
convenient  graphical  form.  Also,  the  stringer's  reactions  have  been 
obtained  for  peak  pressure  intensities  of  PQ  =  2  psi  and  PQ  =  1  psi  for 
different  panel  sizes  and  are  given  in  graphical  form.  An  approximate 
flexural  response  of  the  elliptical  frame  was  formulated,  and  the  maxi¬ 
mum  dynamic  stress  was  determined. 
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II.  ELASTIC  DYNAMIC  ANALYSIS  OF  THE  FUSELAGE 


A.  Determination  of  the  Dynamic  Response  of  a  Panel  Loaded 
Transversely 

The  modes  of  vibration  of  a  plate  with  fixed- fixed  edges  and 
coordinate  system  shown  in  Figure  1  can  be  given  by 


<j>(x,  y)  = 
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Using  energy  approach,  the  frequency  of  harmonic  oscillations  was 
obtained  and  is  given  by 
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where  a  and  b  are  the  dimensions  of  the  skin  panel,  D  the  flexural  rigidity 
of  the  skin  and  mj  is  the  mass  per  unit  of  skin  area. 

In  the  case  of  forced  vibrations,  the  generalized  coordinates 
satisfy  the  differential  equation 


,  t-i 

^mn  +  “miflmn 


(3) 


where  Mmn,  the  generalized  mass  is  given  by 


When  the  distance  between  the  supports  of  the  panel  is  fixed,  the  in-plane 


forces  develop.  Using  virtual  work, the  following  expression  is  obtained 

where  Qmn(t)  is  the  generalized  force  in  the  m,  n  mode,qmn  is  the 
generalized  displacement,  in  the  m,  n  mode,  Wmn  is  the  work  done  by 
the  loading  in  the  m,  n  mode  and  6  stands  for  variation. 

Then 

a/2  b/2 

6Wmn  =  6  J  J  P(x,  y,  t)wmn(x,  y,  t)dxdy 

-a/2  -b/2 


where 

P(x,  y,  t)  is  the  transverse  loading 

wmn(x,  y,  t)  is  the  transverse  dynamic  displacement  of  the  panel 

Ny  is  the  force  developed  at  the  edge  due  to  stretching 
6X  in  the  x- direction 

Ny  is  the  force  developed  at  the  edge  due  to  stretching 
6V  in  the  y-dircction 
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and 


Nx  =.  -f  Eh  ;  Ny 
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but 


Qmn6<W  K  Wmn 

Then 

a/2  -b/2 

Omn^l  ~  J  J  P(x,y,t)4>mn(X,  y)  dxdy 

-a/2  -b/2 


(5) 


In  the  case  under  consideration 

(  2mrrx  \  /  2mTy  \ 

1  +  cos  — • — )  [l  +  cos  -g-  ) 

and  the  Qmn(t)  is  given  by 
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B .  Determination  of  the  Sensitivity  of  Panel  Response  to 
the  Manner  in  Which  the  Panel  is  Loaded 

The  shock  is  travelling  across  the  panel  as  shown  in  Figure  2, 
and  the  loading  is  given  by 

(  < 

P(x,  t)  =  P0\l  -  - 

where 

PQ  =  Peak  pressure 

t  =  Time 

a  -  Dimension  of  panel  in  the  x-direction 
v  =  Velocity  of  shock  propagation 

t^.  =  Duration  of  the  positive  phase 

Using  the  above  P(x,  t),  the  generalized  force  is  given  by 


a/2  +  x 


t+ 


ab 


Ehab 
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q  J  { t) 
min'  ' 


(7) 


The  differential  equation  for  the  generalized  displacement  qmn  is 


„  2  Qmrjt) 

^mn  +  “mnHmn  = 
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where 

Ci^n  is  the  natural  frequency  of  the  panel  in  the  m,  n  mode 
Qmn(t)  is  the  generalized  force  in  the  m,n  mode 
Mmn  is  the  generalized  mass  in  the  m,  n  mode 
For  the  case  under  consideration 
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wmn 


17, 30D 
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where 
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h  =  Thickness  of  the  panel 

E  =  Modulus  of  elasticity  of  the  panel  material 

v  =  Poisson's  ratio 

m 

m  and  n  number  of  half  waves  in  the  x-  and  y-directions  respectively 
m,  =  Mass  per  unit  area  of  panel 


a/2  b/2 

Mmn  =  mi  /  J  'f’rW*’ Y)  dx  =  2-25miab  (8) 

-a/2  -b/2 


The  qmn(tj)  in  terms  of  differences  is  given  by 


^mn(ti) 


^mn^i  +  *>  -  +  qmn^i  -  l) 

(At)2 


where  At 


time  interval 


Substituting  qmn(ti),  in  terms  of  differences,  in  the  differential  equation 
and  solving  for  qmn( +  1)  the  following  difference  equation  is  obtained 
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9mn(ti  ^  1)  -  A{ti)  *  ®mn9mn(ii)  Cmnqmn(ti)  -  qmn(^i  -  1} 


where 
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Considering  the  loading  only  as  a  function  of  time  then 


P(t)  =  P 


and  the  qmn(ti  +  1)  is  given  by 

^rnn^i  f  "  A(tj^)  +  Bmnqmn(tj)  +  Cmnqmn(tj)  -  qlnn(tj  -  1) 


M  tj) 


(At)zP0  , 

2 .  25m  j  \  * 


(9) 


where 


and  the  Bmn  and  Cmn  are  the  same  as  in  the  previous  case, 


The  following  typical  case  has  been  considered 
v  =  1000  ft/ sec 

t+  =1  sec 
F0  =  1  psi 

a  =  20" 

b  =  8" 

h  =  0.03" 

E  =  10.6X10^  (Aluminum) 

v  =0.3 

At  =  0 . 0002  sec 

and 

qm„(0)  =  qmn(0);  (rest  initial  conditions) 

Since  the  1,  1  mode  is  the  prcdom inant  one ,  q^(t)  was  computed 
for  the  two  types  of  loading  and  the  results  obtained  are  given  in  the  table 


below. 


Time  of  qH  max, 

Eoading 

1  max 

sec 

P(x,  t) 

374. 23  X 10'4 

0.0012 

P(t) 

374.  16  X  10"4 

0. 0012 

1 . 00019 


In  the  light  of  the  above  results,  qj  j ,  which  is  a  measure  of 
response  is  essentially  identical  for  both  loading  cases,  not  only  in  terms 


12 


of  magnitude  but  also  in  terms  of  time  dependency.  Thus,  in  the 
remainder  of  the  work,  the  loading  will  be  taken  as  a  function  of  time 
only. 

In  order  to  see  how  the  q.^  (t)  varies  as  the  peak  load  intensity 

varies,  the  response  of  the  previously  considered  panel  was  determined 
for  PQ  of  2  psi  and  4  psi.  The  (t)  were  determined  and  are  given 

in  Figure  3.  The  nonlinearity  is  clearly  shown,  and  the  time  of  maximum 
generalized  displacement  which  is  the  same  as  the  time  at  which  is 
maximum,  as  expected,  becomes  smaller  and  smaller  with  increase  of 
peak  load  intensity. 

C.  Linear  Panel  Response 

In  case  the  stretching  of  the  middle  plane  of  the  panel  is  neglected, 
the  nonlinearity  is  removed  and  the  difference  equation  becomes 

Wli  +  =  A<li>  +  CmnWh)  -  Vn(‘i  -  D 

The  same  panel  has  been  considered  as  before  with  PQ  =  0.1  psi  and 
initial  conditions  corresponding  to  no  motion.  The  maximum  q^(t)  and 
the  time  it  occurred  are  given  in  Figure  3.  By  comparing  the  results  of 
B  and  C  as  given  in  Figure  3,  it  is  clear  that  the  greater  the  PQ  (peak 
load  intensity)  becomes,  the  greater  the  reduction  factor  in  maximum 
displacement  becomes  when  stretching  of  the  panel  is  taken  into  account. 
When  membrane  forces  are  neglected,  the  relation  of  max  and  P0 


is  linear 
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P1  o 1 
o  'Hi  max 

P*  =  q* 
o  11  max 


When  the  membrane  forces  are  included,  the  PQ  vs  <ln(t)inax  is  not 
linear  and  approaches  the  relation 


q' 

11  max 
1 1  max 


D.  Determination  of  Maximum  Direct  and  Bending  Stress 

Utilizing  the  strain- displacement  and  stress-strain  relations, 
the  maximum  direct  stress  in  the  panel  is  given  by 


„  _  2  2 

d  “  ,2  rn  n 


and  the  bending  stress  at  x  =  0,  y  -  y  is  given  by 


o-b  = 


8ir2Eh3 
12(1  -v2) 


m  n  n  qmnH) 


(U) 


(12) 


The  qjj(t)  is  the  predominant  component  and  using  only  this  component 


the  ratio 


is  equal  to 


12(1  -  v‘)qn 


10.1qH 


■=0.3 


For  the  typical  case  considered  in  B,  the  h  =  0. 03"  and  the  max  ^or 


1  psi  was  found  equal  to  about  374.  24  X  10"H  For  this  case: 
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10. 9  X  374  24  X  10’4  10. 9  X  374. 24  X  10-4 


(0. 03) 


27  X  10' 


=  15,100 


Even  if  h  =  0. 05"  and  say  q^  j  max  ^roPs  to  37  X  10”4,  then  the 


^d  _  10  9  X  37  X  10' 

°b  (0 . 0 5 ) 3 


=  322 


E.  Simplification  of  the  Procedure  for  Determining  Panel  Response 
On  the  basis  of  the  above  results,  it  is  clear  that  when  the  in-plane 

forces  are  taken  into  account  the  panel  behaves  primarily  like  a  membrane, 

\ 

which  is  equivalent  to  saying  that  the  bending  stresses  are  insignificant, 
and  the  panel  can  be  considered  as  simply  supported  with  the  distance 
between  its  supports  fixed.  Thus,  for  the  cases  to  be  treated  in  the  future, 
the  panel  will  be  considered  as  a  membrane, and  the  stretching  of  the 
middle  plane  will  be  taken  into  consideration.  The  above  simplification 
'Will  result  Tr T  simplification  of  the  computation  work  without  any  important 
loss  of  accuracy.  | 

F.  Response  of  the  Panel  and  Resulting  Stresses 

1.  Forced  Vibrations  of  a  Membrane  With  Fixed 
Distance  of  Its  Edges 

The  coordinate  system  considered  is  given  in  Figure  4. 

The  dynamic  displacement  of  the  membrane  can  be  represented  by 


w(x,  y,  t)  =  EEqmn(t)<t>mn(x,  y) 
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where  <f>mn(x,  y)  are  the  modes  of  free  vibration  of  the  membrane  and  are 
given  by 


,  y  .  .  nvtTX  niry 

y)  =  91n  — “  8in  -r2- 


.  (13) 


and  <lmnU)  satisfy  the  differential  equation 


2  mn 
^mn  +  wmn<lmn 


M 


mn 


(14) 


where 


CJ. 


mn 


Srr^  /  m2  n2  \ 

\7  +  b1/ 


(15) 


=  Initial  tension  of  the  membrane 


in; 


=  Mass  per  unit  area  of  the  membrane 


Q 


mn 


=  Generalized  force  which  is  given  by 

=  J  J  P(x,  y,  t)  4>mn(x,  y)  dxdy 
0  0 


¥^n(t)/r/  (• 

o  L  o 


^inn^  /*a  (  ^mn^  dx 

-TTJ  dxJ  \-dT~) 

\  o 


dy  (16) 


Eh  3  fh  ( V+mn'f  .  /»b  (  ^mnV 

2b  qmn(tU  J  [  a y  )  dy  J  \  9y  / 

o  LO  o 


dy 


dx 
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ra  rh 

4  I  J  4>mn<x*  V>dxdy 

U  0 


The  Qm  n(t)  and  Mmn  for  the  4,rnn(x*  y)  under  consideration  are  given 


below. 


a  b 

The  /  /  P(x,  y,  t)<t>mn{x,  y)  dxdy 

0  0  P(x,  y,  t)  =  P(t) 


.a  „b 


/t  mtrx  nny  ab 

J  sin  — — —  sin  — —  dxdy  =  P(t)— ^ ^  ( 1  '  cos  m-rr)( 1  -  cos  nrr) 

oo  a  ^nTr 


For  P(t)  =  P, 


Then 


.a  b 


ra  ru  p°ah  (y  t  \ 

J  J  F(t)4>nin(x.  y)dxdy  -  - — j  ( 1  -  cos  mir)(  1  -  cos  rnr)  \  1  -  —  J 

n  n  i  ■  t 


0  0 


a  .b 


but  for  m  and/or  n  even,  the  I  J  p(^mn(x-  l)dxdy  r  0 


0  0 


.a  .b 


'■  /  /  P(t)«j>mn(x,  y)  dxdy 


4abP 


0  0 


rO-77) 


m,  n  =  odd 


Proceeding  in  a  similar  way,  the  other  integrals  ol  (16)  arc  computed, 


and  finally  the  generalized  force  is  given  by: 
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4abP 


m,  n  =  odd  mniT 


:“('-rr)-^EhabK  +  S]q-(t) 


0. 406abP 


m,  n  =  odd 


—  ^  1  -  j-j  -  4.  55  Ehab  £ 


—  +  —  q3  (t)  (18) 

a 4  b4J  mn 


The  generalized  mass  is  given  by 


.a 


a  b 


n  =  mi  /  /  4n(x-y)dxdy  =  mi/  J  9 


0  0 


0  ,  0 


2  mirx  ?  nftx 
in  — —  sin  ~g~  dxdy 


m  .ab 
i 


=  0.25  m^ab 


For  S  (Initial  Tension  of  the  Membrane)  =  0 

Equation  (15)  yields  Wfnn  =  0>  ar>d  the  Equation  (14)  becomes 


^mn(l)  = 


In  terms  of  differences 


'Imn^i  "  1)  ’  ^mn^i)  +  ^mn^i  +  ^ 
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and  Equation  (20)  reduces  to 


^mn^i  "  mn^  *i  +  ^ 


0.  406abP 


abP^  /  t  \  \  r~4  4“ 

^(..^)-4.55Ehab[^+^]q^nUi, 


0.  25m^ab 


‘Imn^i  +  1)  ~  ^mn^i)  +  ^mn^i)  +  Vn^mn^i)  "  ’Imn^i  ”  *) 


where 


Amn<‘i>  = 


1. 624(  AtrP. 


mjmn 


<-*> 


18.  20Eh(  At) 


2  r  4  4*1 

m  n 

"  b3. 


2.  Determination  of  the  Maximum  Stress  in  the  Panel 


The  in-plane  forces  Nx  and  Ny,  as  shown  in  Figure  4,  are 


given  by 


xmn 


6x  Eh  (• *  (' 

—  Eh  =  —  1/2/  ^ 


Using  4>mn(x,  y),  the  expression  given  by  Equation  ( 1 3),  the  (22)  becomes 


( mTr\  .  2  n-rry  ?  i 

{—)  sm  -r  <,™|t) 
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'and  the  stress  due  to  Nxmn  is  given  by 


N 

^xmn 

=  -( 

rmr 

2  \ 

)  aZ 

max  ^ 

max  4  ' 

^  a 

)  qmn^  ^ 

The  same  way 


ymn 


max 


ymn 

h 


max 


4(b)  q™n(t) 


The  maximum  stress  takes  place  in  the  shorter  direction  and  it  is 
equal  to 


<r  =  <r 

mn  max  mn 


/  \2  « 

2 

II 

•MW 

)  qmn^)  =2.46e(- 

1  x  =  a/2 

max 

max 

=  |rSZ.46E(f)‘:q^n,t)j 


Using  Equation  (21),  the  response  of  the  panels  was  determined,  and, 

'  \ 

employing  Equation  (23),  the  maximum  sti^sses  were  obtained  and  are 
given  in  nondimensional  form  in  Figvire  5.'  * 

G.  Determination  of  the  Load  on  the  Stringer 

\ 

According  to  Figure  4,  the  reaction  at  ji  =  0  is 

I 


R  =  N 

x=  0 


Tmn 

xmn  9v  qmn^ 

x  =  0 


359 
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but 


N. 


Eh 
xmn  4 


[srf"Js?4>  n(« 


and 


3<t>m 

_  ( 

1  mTTx  niry 

( mTT  ^ 

|  niry 

Ox 

V  a  ) 

x  -  0 

1  cos  sin  , 

a  b 

-  V  a  i 

x  -  0 

!  sin  , 
b 

or 


R 


x=  0 


\  .  3  niry  3  ,  , 

sin  b  ^mn^ 


=  SS7.74Eh 


(mr  3  niry  3  ,  , 

t)  6in  b  wl> 


(24) 


The  same  way  the  R|  ,  which  is  the  lotiding  on  the  stringer,  is 

|  y  =  0 

obtained  and  is  given  by 


R 


EE  7. 74Eh 


y-  0 


3  rmrx  3 
sin  ——  qmn(t) 


(25) 


The 


=  S27.74Eh(f)  q^n(t) 


y  -  0  max 

In  order  to  obtain  the  relative  contribution  of  the  to  the  R|  „  , 

^mn  y  =  0  max. 

a  typical  panel  was  considered  (a  =  20",  b  =  8";  h  =  0.03"),  and  subjected 
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to  a  linearly  decaying  dynamic  load  of  peak  intensity  equal  to  Pq  =  1  psi, 
and  the  qmn  were  determined  when  reached  its  maximum  value. 
Their  magnitudes  are 

qn  =  0.1619 
q13  -  0.0005 

q3l  =  0.  0266 

q33  =  0.0183 

Their  contribution  to  the  R|  „  is 

|  y  =  0  max 


<?. 

=  4,250  X  10‘6  # 

(3)  3qf  3 

*  0 

3 

q13 

=  18.  8  X  10'6 

(3)3q^ 

56  X  10'6 

4,324.8  X 10"6 

It  '.s  clearly  shown  that  the  contribution  of  the  j  alone  is  98.  3%  of  the 
total  magnitude.  On  the  above  basis,  the  loading  of  the  stringer  will  be 
approximated  by  the  contribution  of  the  q^(t)  alone,  and  thus  will  be 

given  by 


R 


y=0 


7 . 7  4Eh  -i  ux  7. 
—-y—  sin  —  qf1<t) 


(26) 


3  ^ 

In  order  to  get  an  idea  how  the  j(t)  varies  with  time,  the  q^^t)  has  been 
plotted,  for  the  typical  panel  considered  previously,  and  is  given  in  Figure  6. 
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Generalized  force 
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The  generalized  mass  is  given  by  Mi  =  mi J  <jn  (x)  dx 

0 


or 


Mi  =  mly. 

0 


if  [  Mcos  ki 


x  -  ch  kix)  +  (sin  kiX  -  sh  kix)]  dx 


After  performing  the  integration,  the  Mi  can  be  written  as 


Mi  =  m^  o-i 


(32) 


where 

a-  =  — - —  [  \2(  2k;.(!  +  sin  k:i  cos  kxl  +  sh  k;i  ch  k^i  -  2  shk;^  cos  k^f  -  2chk;i  sink^i)  + 
1  2k;£  x  i  i  i  i 

+  (  -  sin  k^f  cos  k^f  +  sh  k^  ch  kji  -  2  ch  kjf  sin  kjf  +  2  sh  kjf  cos  kji  )  + 

+  2Msin2  kjl  -  2  shk^  sink^  +  sh2  kti  )]  (33) 


Since  the  generalized  force,  Equation  (30),  is  given  in  tabular  form,  then 
the  solution  of  Equation  (29)  has  to  be  obtained  by  numerical  or  graphical 
means. 

I.  Stringer  Response 

Due  to  continuity  of  the  stringer  at  the  frame,  the  ends  of  the 
stringers  will  be  considered  fixed.  The  transverse  displacement  of  the 
stringers  is  much  greater  than  that  of  the  frame,  and,  for  this  reason, 


the  motion  of  the  fixed  ends  of  the  stringer  will  be  neglected  in  obtaining 


I 


7' 

1* 


)i 

L 


r 


i; 
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the  stringer  response.  The  dynamic  displacement  of  the  stringer  is 
given  by 


*(x,t)  =  S^xjq^t) 
i 


where 


<jq(x)  =  (X^cos  k^x  -  ch  kj,x)  +  (sin  k^x  -  sh  k^x)] 


The  generalized  coordinates  q^(t)  satisfy  the  differential  equation  and  the 
initial  conditions  given  below 

2  °i(t) 

qi(t)  +  ufqitt)  =  ,  qi<0)  =  qi(0)  =  0 


(34) 


qi  i(tM  ( qi  i(t) 

The  Oj.{t)  =  7 . 7 4Ehapi  ^ ^ — )  ,  and  the 


is  given  in 


graphical  form  from  the  skin  response  in  Figures  8,  9,  and  10.  Dimen¬ 
sional  analysis  shows  that  the  time  scales  as  b  (i.e.  ,  kt  =  k^).  and 
Figures  8,  9.  and  10  can  be  used  for  the  same  a^h/b^  for  the  same  peak 
pressure,  PG,  and  density,  p,  of  the  skin  material.  Equation  (34)  has  to 
be  integrated  numerically  or  graphically. 

In  integrating  Equation  (34),  a  numerical  procedure  will  be  used 
and  its  formulation  is  established  below.  The  qjjt)  in  terms  of  differences 
is  given  by 


9i(ti  -  1)  -  2qi(ti)  +  qi(ti  +  1) 
(At)2 


34 


then  Equation  (34)  becomes 


qi(ti  -  1)  -  2qi(ti)  +  qi<4  4  1)  =  - 


(*tf  4-0^)  +  - jjp — 


or 


Qi(tj)(At)2  ,  , 

qi(ti  +  1)  =  - -  +  l2*  <At>‘ "i  1  -  1) 


or 


qi(t*  +  1)  =  A^tj)  +  Bjq^ti)  -  q^tj  -  1) 


(35) 


where 


Ai(ti) 


3 

7 . 7  4Eh  (  qli(ti)  )  7  7  74Eh  1 

=  —)  (Atr=“^~ pi' 


Bj  =  (2  -  (At)**) 


wi 

E 

h 

m; 


q,,(t) 


2  /eT  ! 

=  (k-  ) — ■  =  Frequency  of  the  stringer 

=  Modulus  of  elasticity  of  the  stringer's  material  =  10.  6X  103  psi 
=  Thickness  of  skin 
= Mas s  per  unit  length  of  stringer 
=  LenRth  of  stringer 
=  Width  of  panel 

=  Generalized  coordinate  of  the  skin’s  response  in  the 
1  ,  1  tli  mode 


(At) 

=  Time  interval 

kj  a 

=  4.730; 

al  = 

1.244  ; 

Pj  =  0.43682  ; 

X  =  -1.0175 

k2a 

=  7.853; 

a->  = 

L* 

1.00088  ; 

P2  =  0.001212; 

X  =  -1. 0000 

k^a 

=  10.993; 

P 

II 

1. 000000; 

p3  \- -0. 110407  ; 

X3  =  -0.99997 

The  moment  on  the  stringer  is  given  by 

v 

\ 

M  =  -EI®"(x,  t)  =  -  EIEq.(t)<t>’.'  (x)  \ 

1  1  1  '  \ 

-  -  El  £  qt(  t)  kf  [\j{-  cos  kjx  -  chkjx)  +  (-sin  kjx  -  sh  kjx)] 

________  \ 

The  maximum  moment  was  found  to  take  place  at  x  =  0  and  thus 


Mmax  =  M|x=0 


El  Eqi(t)kf(-2X.)  =  2EI  Eq^OkfX. 


a  i 


The  shear  is  given  by 


V  =  -  El®' 


El  EqjO)  4>im  (*) 
i 


=  -  El  Eqj(t)  kf  l  \j( sin  kjx  -  sh  kjx)  +  (-  cos  kjx  -  ch  kjx)] 


The  stringer's  reaction  is  equal  to  V 


x  =  0 
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V 


x  —  0 


-  El  Eq^t)  k?  (-2)  =  2EISqi(t)kf 
l  i 


#  • 


V 


x  -  0 


’El  3 

“  Sq^M^a)3 

a-3  i 


(37) 


Stringer  cross  sections  that  are  considered  to  be  typical  are 
shown  in  FigvrrrTTTr^'The  (a)  has  been  used  in  computing  the  stringer 
response  using  Equation  (35).  The  effective  mass  per  unit  length  of 
stringer  is  not  easily  determined  because  of  the  unknown  interaction  of 
the  skin.  For  this  reason,  first,  the  response  of  the  stringers  was 
determined  without  any  interaction  of  the  skin  (i.  e.  ,  m^  =  mass  per 
unit  length  of  the  stringer),  and  the  results  are  given  in  Figures  12  and  13. 
Also,  computations  were  made  with  full  interaction  of  the  skin,  as  shown 
in  Figure  14,  which  resulted  in  increasing  the  effective  mass  per  unit 
length  of  stringer  by  a  factor  of  about  three.  The  response  of  the  stringer 
was  determined,  and  the  obtained  moments  and  stress  were  found  about 
30%  smaller  than  without  any  skin  interaction. 

Because  of  the  fact  that  an  interaction  of  skin  and  stringer  exists 

I 

and  is  between  complete  and  zero  interaction,  it  is  believed  that  such  an 

interaction  will  reduce  the  response  given  in  Figures  12  and  1  3  by  15%. 

I 

Thus,  if  the  values  given  in  Figures  12  and  13  are  modified  by  -15%,  they 
are  going  to  represent  the  true  response  of  the  stringer  to  within  5%, 
and  this  is  considered  satisfactory  for  the  problem  at  hand.  Such  modified 
curves  lor  direct  use  in  design  are  given  in  "Figures  15  and  16. 
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FIGURE  i  1 .  TYPICAL  CROSS  SECTIONS  OF  STRINGERS 
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FIGURE  I4-.  SCHEMATIC  REPRESENTATION  OF  FULL  SKIN  INTERACTION 


b  FOR  P  =  2  PSI,  AND  LOADING  OF  ONE  PANEL 
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Loading  curves  as  well  as  modified  curves  for  direct  use  in  design 
have  been  obtained  for  peak  pressure  intensity  PQ  =  1  psi  and  are  given 

in  Figures  17  through  21. 

Based  on  the  discussion  that  follows,  the  results  presented  in 
Figures  15,  16,  20,  and  21  can  be  used  for  other  stiffeners  as  well. 

For  a  stringer  with  the  same  I/rm  as  the  one  used  herein,  the 
frequency  will  be  the  same,  and,  for  an  increase  of  rm  by  50%  and  even 
100%,  the  moment  and  shear  will  be  proportional  to  I  and  not  off  by  more 
than  -10%.  The  bending  stresses  are  the  critical  ones  in  a  stringer;  thus, 
for  the  same  I/m^.  as  the  one  used  herein  and  an  increase  of  m^  by  50%  or 
100%,  the  bending  moments  will  be  proportional  to  I  and  not  off  by  more 
than  -10%.  In  this  way,  the  results  of  Figures  15,  16,  20,  and  21  can  be 
used  for  design  purposes  when  determination  of  maximum  stress  in  the 
stringer  with  less  .than  ±10%  off  can  be  tolerated. 
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J.  Formulation  of  the  Approximate  Flexural  Dynamic  Response 

of  an  Elliptical  Ring 

The  differential  equations  of  motion  of  an  elliptical  ring  of  constant 
cross  section  are  linear  of  the  fourth  order  with  highly  complex  variable 
coefficients.  Their  solution  in  terms  of  known  functions  -is  inordinately 
difficult  and  is  considered  t.o  be  beyond  the  scope  of  this  investigation. 

Thus,  the  true  modes  of  vibration  have  not  been  determined. 

Because  of  that,  an  approximate  solution  is  sought  on  the  basis  of 
an  assumed  mode  of  fundamental  vibration.  A  typical  elliptical  frame  of  a 
fuselage  has  been  provided  and  is  given  in  Figure  22.  It  :.s  true  that  the 
cross  section  of  the  frame  is  not  constant,  and  also  that  about  seventy- 
five  percent  of  its  length  has  a  constant  cross  section.  Also,  the  frame 
under  periodic  loading  along  the  circumference  will  deflect  more  in  the 
direction  of  the  minor  axis  than  in  the  direction  of  the  major  axis  of  the 
ellipse.  This  difference  in  the  magnitude  of  the  above  mentioned  displace¬ 
ments  becomes  greater  for  greater  ellipticities .  The  increase  in  rigidity 
at  the  ends  of  the  major  axis  will  result  in  a  decrease  of  the  displacements 
in  that  portion  of  the  frame,  and.  at  the  same  time,  will  decrease,  in  a 
smaller  degree,  the  displacements  of  the  frame  where  the  cross  section 
remains  constant.  It  is  believed  that  the  ratio  of  the  areas  under  the 
static  deflection  curve  for  a  frame  with  constant  cross  section  and  the 
same  ratio  of  the  areas  under  the  deflection  curve  for  a  frame  with  increased 
rigidity  at  the  ends  of  the  major  axis,  under  the  same  loading,  and  whose 
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FIGURE  ’.2.  TYPICAL  FRAME 
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seventy-five  percent  of  the  length  has  a  constant  cross  section,  differ 


*  . .. 


but  little.  If  the  dynamic  response  of  the  constant  section  and  variable 
section  frame  is  considered  on  the  basis  of  their  statit  deflection  curve 
mode,  the  generalized  force  {i.  e.  ,  load  times  mode  displacement)  for 
both  frames  will  be  essentially  the  same  for  both  frames.  If  the  generalized 
mass  is  considered  (i.  e.  ,  mass  times  mode  displacement  squared)  Jt  will 
differ  very  little  for  the  two  frames, because  the  increase  in  mass  takes 
place  only  on  twenty-five  percent  of  the  frame  length,  and  because  the 
displacement  in  that  portion  of  the  frame  is  much  smaller  than  in  the  flat 
portion  of  the  ellipse.  In  addition,  the  frequency  of  the  constant  section 
and  variable  section  frames  will  be  different  by  a  small  amount  which  will 
cause  only  a  slight  shift  in  the  time  of  maximum  response,  leaving  the 
magnitude  of  the  generalized  displacement  slightly  affected.  Based  on 


the  above  discussion,  it  is  believed  that  under  the  loading  at  hand  the  frame 

with  localized  increase  in  rigidity  and  one  with  coinstant  cross  section, 

\ 

equal  to  the  one  existing  on  the  seventy-five  percertt  of  the  frame  with 

\ 

variable  cross  section,  will  have  practica’ly  the  same  deflected  shape 

% 

and  dynamic  response. 

It  is  known_that  any  deflection  curve  can  be  expressed  in  tqrms  of 


the  natural  modes.  In  the  absence  of  the  true  modes,  it  is  believed  that 


the  best  approximate  shape  to  be  used  as  a  mode  would  be  the  static 


deflection  curve. 
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The  loading  of  the  frame  is  concentrated,  coming  from  the  stringers' 

\ 

reactions,  and  distributed,  coming  from  the  skin.  Due  to  relatively  close 
spacing  of  the  stringers,  the  static  deflection  curve  under  this  loading 
will  be  close  to  the  static  deflection  curve  of  the  constant  section  elliptical 
ring  under  uniform^loading .  Such  a  deflection  curve  is  not  difficult  to 

obtain,  and  for  an  approximate  treatment  it  is  considered  satisfactory. 

\ 

As  a  result  of  the  above  reasoning,  the  approximate  flexural 
dynamic  response  of  the  elliptical  ring  at  hand  will  be  based  on  the  static 
deflection  curve  mode  under  uniform  loading.  Such  a  mode  is  shown  in 
Figure  23. 

If  it  is  denoted  by  <f>(s),  where  s  is  the  measure  of  length  of  the 
ellipse,  theri"the""dynamic  displacement  is  given  by 

\ 

«(s,t)  =  <j>{s)  q(t)  (38) 

where  the  q(t)  satisfies  the  differential  equation  and  initial  conditions 

q(t)+co2q(t)  =  ;  q(0)  =  q(0)  =  0  (39) 

where  w  =  the  frequency  and  is  given  by 


1  0  to  tt/2  due  to  symmetry 
of  4>(s)) 


^&WeC0PV 


duced  F'°m 


(40) 
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(41) 


The  M  =  generalized  mass,  which  for  constant  cross  section  is  given 


o  r 


\ 


62 


=  -  Ec 


max 


az4>(s) 

as2 


<WxW 


(44) 


max 


K.  Determination  of  the  Dynamic  Response  of  a  Typical  Elliptical 

Frame 

The  typical  frame  considered,  with  the  attached  stringers  and  skin, 
is  given  in  Figure  30.  Considering  the  loading  coming  from  one  side  of 
the  frame,  du^^io^intensity  of  2  psi  and  with  positive  phase  duration 
of  1  sec,  the  maximum  bending  stress  was  detcrjprnned  to  be  equal  to 
crj  -  93,  200  psi 

L.  Comparison  of  the  Order  of  Magnitude  of  the  Static  and  Dynamic 

Stresses  and  Additional  Static  Deflection  Curves 

The  determination  of  the  maximum  bending  moment  in  an  elliptical 
ring  under  uniform  static  pressure  is  given  by  S.  Timoshenko.*  This 
moment,  for  the  geometry  considered  in  the  previous  case,  is  given  by 
0.870  Pb2,  where  P  is  the  uniform  load  intensity  on  the  frame,  and  b  is 
half  of  the  minor  axis  of  the  ellipse. 

If  P  =  (2)(9)  =  18  pli  (i.e.  ,  half  of  the  panel  load  per  length  of 
frame),  then  the  maximum  static  bending  stress  would  be 


=  0.  870(18)(16)2 


1.  505 
0.147 


41 1 000  psi 


*S.  Timoshenko,  "Strength  of  Materials ,  "  Part  II,  Second  Edition,  p.  90. 
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FIGURE  30.  TYPICAL  FUSELAGE  SECTION 
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The  same  frame  was  considered  under  uniform  dynamic  pressure  of  peak 

intensity  PQ  =  18  pli,  and  the  maximum  dynamic  stress  obtained  was 

found  equal  to  73,000  psi.  Thus,  the  ratio  of  the  maximum  dynamic 

to  maximum  static  pressure  for  peak  intensity  equal  to  the  static  load- 

73,  000 

intensity  is  equal  to  ^"'"qqq  =  1-78.  It  should  be  observed  that  the 

frame  under  consideration  is  not  loaded  uniformly,  and  that  most  of  tho 
loading  is  a  concentrated  loading  coming  from  thd  stringers'  reactions. 
For  such  a  loading,  the  would  be  greater  than  7j3,000  psi,  and,  thus, 
compared  with  the  dynamic  bending  stress  obtained  in  the  previous 
section  (i.e.  ,  cr^  =  93,200  psi)  results  in  a  check  on  the  order  of 
magnitude  of  the  obtained  maximum  dynamic  bending  stress. 

Static  deflection  curves  have  been  obtained  for  a/b  equal  to  2.  5 
and  1. 666.  They  are  given  in  Figures  31  and  32  and  can  be  used  for  any 


a  and  h  of  the  same  ratio. 


rwimss 
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M .  Sample  Calculations 

In  this  section,  the  response*  of  the  skin,  stringer,  and  frame 
shown  in  Figure  30  will  be  considered  in  detail  for  peak  pressure  intensity 
of  PD  -  2  psi,  and  positive  time  duration  t  +  =  1  sec. 

^  1.  — i S’WTf"fresponse 

The  generalized  coordinate  of  the  skin  response  is  governed 
by  Equation  (21)  given  below 

^mn^i  +  1)  =  +  ^mn^i)  +  ^mn^mn^i)  "  ^mn(li  "  M 

where 


^mn^i) 


1. 624(At)2P( 
mjmn 


(■-*) 


‘•'mil 


18,  20Eh(At)Z 


For 

E  -  10.  6  X  106  psi 

h  =0.03  in. 

PQ  =  2  psi 

At  =  0 . 000 1  sec 

a  =18  in. 

b  =4  in. 
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=  7 . 42  X  10'6  lb  in’1  sec2 


=  -  30. 5435 


-  2,468.068 


=  -  36.4877 


=  -  2,  474.  0196 


The  Amn(t  j)  have  been  computed  and  are  given  in'tabular  form  in 
Tables  1,  2,  3,  and  4  where  the  numerical  integration  of  Equation  (21) 

is  performed  for  the  determination  of  the  generalized  coordinates  qmn(t) 

\ 

The  maximum  stress  in  the  skin  is  given  by  Equation  (23) 


ffraax  =  |rS2.46E(i)2 


max 


or,  usingthe  results  obtained  in  Tables  1,  2,  3,  and  4,  the 


max 


(10. 6  X  106)  f  ,  ?  ?  ?  , 

- 5 -  1  (0. 08376)2  +  (3r  (0.  00523)2  +  (0.  03865)2 

(4)  l  t 


+  (3)2  (0.  009306)2j-  =  15,600 


2 .  Stringer  Response 

The  generalized  coordinate  of  the  stringer  response  is  given 


by  Equation  (35) 


+  l)  =  Ai<li)  +  B^ijdj)  -  qjdj  -1) 


TABLE  1.  DETERMINATION  OF  qn(t) 
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Ai«i) 


rrijQj  \  b  / 


Bj  =  [2  -  «f  (At)2] 


The  A^(t)  are  practically  proportional  to  Pi  and  thus  the  A£(t)  is  negligible 
compared  to  A^t).  The  A^t)  and  Az(t),  as  well  as  the  Bi  and  are 

computed  and  given  in  Tables  5  and  6  with  the  nurherical  integration  of 

\ 

Equation  (30)  for  the  determination  of  q j ( t)  and  q3< t) .  The  ^j(t)  and  q${t) 
are  given  in  Tables  5  and  6,  and' it  is  clearly  shown  that  the  contribution 
of  qi(t)  in  determining  the  dynamic  response  of  the  stringer  can  be 

neglected,  and,  by  doing  that,  the  response  obtained  is  slightly  on  the 

\ 

\ 

conservative  side. 

Thus,  the  maximum  bending  moment  and  shear  at  x  =  0  as 
given  by  Equations  (36)  and  (37)  reduce  to 


^  'll  max't>'kla>2M 


» and  , 


max 
x=  0 


2EI 

T  ** 


1  max 


(t)(kja) 


TABLE  5.  DETERMINATION  OF  q  (t) 
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TABLE  6. 


DETERMINATION  OF  q3  (t) 


t-  SGC 

(1) 

A3Ui) 

(2) 

q3(ti) 

(3) 

1. 6207  (2) 

\  (4) 

qi(ti  -  i) 

(5)  =  (1)  +  (3)  -  (4) 

q3(  t  i  +  i) 

0 

0  X  10"6 

0  X  10“6 

o  x  io-6 

oix  10-6 

•  0  X  10-6 

0. 00020 

0.  2544 

0 

0 

1i> 

0  .  259  4 

0.00022 

0.  44 

0. 2594 

0. 4204 

0 

0.8604 

0. 00024 

0.  64 

0.  8604 

1.  394 

0. 25S4 

\  1.7751 

0. 00026 

0.  84 

1. 7751 

2.  877 

0. 8604 

'  2.857 

0. 00028 

1. 04 

2.857 

4.  6303 

1 . 7751 

3.  895 

0. 00030 

1 . 206 

3.  895 

6.  3126 

2.  857 

4.  6616 

0.00032 

2.000 

4. 6616 

7.555'  ' 

3.  895 

5. 6600 

0. 00034 

3.2 _ 

5.  6600 

9-  1732 

4.  6616 

7  7116 

0.00036 

4.  6 

7. 7116 

12. 4982 

5.  6600 

11. 4381 

0. 00038 

6.  8 

11. 4381 

18. 5377 

\  7.7116 

17. 6261 

0. 00040 

9.  16 

17.6261 

28. 5666 

'll  4381 

26. 289 

0. 00042 

12.  00 

26. 289 

42. 6066 

17. 6261 

36.9805 

0. 00044 

14.  80 

36. 9805 

59. 934 

26. 289 

48. 445 

0. 00046 

18.  40 

48. 445 

78. 5148 

36.  9805 

59. 934 

0.00048 

22.  40 

59. 934 

97. 135 

48.  445 

71.090 

0. 00050 

27.  25 

71. 090 

115.  216 

59.934 

82. 532 

0. 00052 

32.  80 

82. 532 

133.  760 

71. 090 

95.  470 

0.  000  54 

38.  4 

95.  *70 

1  54.  728 

82  532 

110. 596 

0.00056 

43.  2 

110.  596 

179. 243 

95.  470 

126.  973 

0. 00058 

49.  2 

126. 973 

205. 785 

1 10.  596 

144.  389 

0. 00060 

144. 389 

234  011 

126. 973 

161 . 248 

0. 00062 

58.  40 

161 . 248 

261. 335 

144.  389 

l*71.  345 

0.00064 

62.  4 

175.  345 

284. 1816 

1161 . 248 

185.  334 

0.00066 

64.  8 

185  334 

300. 371 

175.  345 

189. 826 

0.00068 

67. 2 

189. 826 

307. 691 

185.  334 

139. 517 

0.00070 

67.  6 

189.  517 

307.  1  50 

189. 826 

184.  924 
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P'or  the  case  under  consideration,  the 


M 


2(10.  6  X  106)  I  (0. 0244)(4.73)Z(rl.  0175) 


max 


(18)' 


f 


—  =  364  X  102  I 


Me 


max 


—  =  364  X  10^(0.614)  =  22,  300  psi 


and 


(from  the  loading  of  one  panel) 


x=  0 
max 


2(10.  6  X  106)  ? 

=  — - x - -  (0. 0095)(0.  0244)(4.  73)3  =  89-6  lb 


(18) 


(from  the  loading  of  one  panel) 

3.  Frame  Response 

For  the  frame  under  consideration,  t^e  mode  shape  to  be 

3 

used  has  been  determined  and  given  in  Figure  33.  pbtaining  (qjj/a) 

V. 

from  Figure  24,  qj(t)  from  Table  5,  using  the  mode  given  in  Figure  33, 

\ 

and  using  numerical  integration,  thew,  Q(t)  and  M  given  by  Equations  (40), 
(41)  and  (42)  have  been  computed,  and  then  Equation  (43)  was  integrated 
to  determine  the  generalized  coordinate  q(t) 

The  numerical  integration  of  Equation  (43)  is  given  in 

\ 

tabular  form  in  Table  7.  Then  the  maximum  bending  stress  is  given  by 


the  Equation  (44): 
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TABLE  7.  DETERMINATION  OF  q(t) 
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The  - j—  is  maximum  at  s  =  0.  The  origin  of  the  curve  on  Figure  33 

ds* 

lias  been  replotted  in  larger  scale  and  is  given  in  Figure  34.  Numerical 
differentiation  yields 


Then 


9Z4)(o)  24»(Ab)  -  2<fr(o)  2(.Q.  080)  -2(-0. 0825) 


9  s' 


(As)- 


(0.  65)2 


0.01185 


<r_=10.6  X  10°(1.505)(0.01185)(0.4925)  =  93,  200  psi 

m  aX 

(from  the  loading  of  one  side  of  th^  frame) 
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